A study of heuristic guesses for adiabatic quantum computation 



O 
(N 

(N 






> 

m 

p 

o 

oo 
O 



X 



Alejandro Pcrdomo,^ Salvador E. Venegas-Andraca,^' ^ and Alan Aspuru-Guzik-'^ 

'Department of Chemistry and Chemical Biology, 

Harvard University, 12 Oxford Street, 02138, Cambridge, MA 

'^ Quantum Information Processing Croup. Tecnologico de Monterrey Campus Estado de 

Mexico. Carretera Lago Cpe. Km 3.5, Atizapdn de Zaragoza, Edo. Mexico, Mexico 

Adiabatic quantum computation (AQC) is a universal model for quantum computation which 
seeks to transform the initial ground state of a quantum system into a final ground state encoding 
the answer to a computational problem. AQC initial Hamiltonians conventionally have a uniform 
superposition as ground state. We diverge from this practice by introducing a simple form of 
heuristics: the ability to start the quantum evolution with a state which is a guess to the solution 
of the problem. With this goal in mind, we explain the viability of this approach and the needed 
modifications to the conventional AQC (CAQC) algorithm. By performing a numerical study on 
hard-to-satisfy 6 and 7 bit random instances of the satisfiability problem (3-SAT), we show how 
this heuristic approach is possible and we identify that the performance of the particular algorithm 
proposed is largely determined by the Hamming distance of the chosen initial guess state with 
respect to the solution. Besides the possibility of introducing educated guesses as initial states, the 
new strategy allows for the possibility of restarting a failed adiabatic process from the measured 
excited state as opposed to restarting from the full superposition of states as in CAQC. The outcome 
of the measurement can be used as a more refined guess state to restart the adiabatic evolution. 
This concatenated restart process is another heuristic that the CAQC strategy cannot capture. 

PACS numbers: 03.67.Lx, 03.67.Ac, 03.65.-w 



I. INTRODUCTION 

Adiabatic quantum computation (AQC) [l| is a 
promising paradigm of quantum computation because of 
its robustness [2|, |3| , and its intuitive mapping from NP- 
complete and NP-hard problems to potentially realizable 
Hamiltonians [l|, |j-|6| . Adiabatic quantum computing is 
attractive because relevant optimization problems such 
as lattice models for protein folding can be readily for- 
mulated 0- 

AQC algorithms involve the specification of a time- 
dependent Hamiltonian, 



H{t) = hi{t) + hdriving{t) + h f {t) , 



(1) 



Since the original proposal by Farhi et al ]l[, a sig- 
nificant amount of progress has been made towards the 
design of final Hamiltonians for different computation- 
ally intractable problems such as NP-complete problems. 
For example, in this paper we provide a detailed descrip- 
tion of the construction of final Hamiltonian for the 3- 
SAT problem. In previous work, we described the re- 
spective construction of the final Hamiltonian for an NP- 
hard problem of interest in biology, the protein folding 
problem, which consists of finding the minimum energy 
configuration of a chain of interacting amino acids in a 
lattice [Tj. Several choices for hdriving have been sug- 
gested. Farhi et al [l[ proposed what we call in this work 
the conventional form used for hdriving{t)'- a linear ramp, 



This Hamiltonian has three important functions: (1) 
The initial Hamiltonian, Hi = H{0), encodes a ground 
state that is easy to prepare and that is used as the ini- 
tial state for the quantum evolution. (2) The driving 
Hamiltonian, hdriving (t), is responsible for mediating the 
transformation of the initial ground state to any of other 
state. (3) The final Hamiltonian, Hf = H{t), is prob- 
lem dependent and its ground state encodes the solu- 
tion, \ijJsoiution) , to the computational problem. In the 
ideal case of a process being fully adiabatic, evolution 
under H(t) will keep the quantum state, \Tp(t)), in the 
ground state of H{t) throughout < t < t. If this 
condition is met, the final state a.t t — t should coin- 
cide with the ground state of the final Hamiltonian, Hf , 
i.e., |'0(t)) = \ipsoiution) , if the process is adiabatic. The 
measurement at t = t will provide the solution to the 
computational problem. 



'^drivingy^) — yt Z/TjTldr 



(2) 



In Farhi et aVs scheme the term hi{t) can be defined 
as iJj = hdrvving{t = 0) = H driving, namely, the term 
Hdriving scrvcs as the initial Hamiltonian Hi when its in- 
tensity is completely turned on at t = 0, and drives the 
quantum evolution until it is completely turned off at 
t — T. The final Hamiltonian is turned on with the func- 
tional form hf{t) = {t/T)Hf. Scheduling the adiabatic 
evolution with this linear interpolation is not compulsory, 
thus different proposals have been studied such as the use 
of non-linear time-profile for auxiliary Hamiltonians [8| 
and optimal geodesic paths Q- 

The question of whether AQC, or in general quan- 
tum computers, can be used for efficiently solving NP- 
complete problems is a difficult open question. Lessons 
that have been learned [1(| include that a poor choice 



of the initial Hamiltonian such as the one-dimensional 
projector as selected in Refs. jll[ and [l3, will lead in- 
efficient AQC algorithms. Therefore, it is important to 
consider different strategies which might allow an escape 
from bottlenecks or trap states which might limit the use 
of AQC [H. 

When attempting to tackle combinatorial and opti- 
mization problems with classical computers, a common 
approach to cope with intractability and NP-complete 
problems [ij, |l^ , is to employ heuristic algorithms as an 
alternative to exhaustive search which scales exponen- 
tiall y w ith system size. In quantum computation, Hogg 
[la . Il7| introduced heuristic techniques in quantum al- 
gorithms. He showed that using information about the 
structure of the problem as a heuristic guide can be used 
to enhance the performance of quantum search compared 
to the scheme proposed by Grover ^^. Hogg's proposal 
was suggested for the gate model for quantum computa- 
tion, but to our knowledge, there are no studies of heuris- 
tics in AQC. The purpose of this paper is to examine the 
following questions: How can we incorporate heuristics 
in AQC? Is there any advantage by doing so? What are 
the modifications to CAQC proposal from an algorithmic 
and experimental point of view? 

a. Initial state selection. The simplest form of 
heuristics we could think of is to start the quantum evo- 
lution from a quantum state which is a guess to the solu- 
tion, but this possibility is not available in any of the pro- 
posals for AQC. Physical intuition as well as constraints 
within the problem can be used to make an educated 
guess. To illustrate our idea, let's use a lattice model 
for protein folding 7] as an example. For this model, an 
educated guess for the initial state would be to choose a 
bit string which encodes an initial position for the amino 
acids in the spatial lattice such that no two amino acids 
are on top of each other, and that they are connected ac- 
cording to the sequence defining the protein to be folded. 
Conversely, CAQC would begin the computation with a 
quantum superposition of all possible states of the com- 
putational basis. This choice of initial state contains ab- 
surd configurations like all amino acids on top of each 
other, or assignments referring to configurations of amino 
acids which are fully disconnected or not properly linked 
according to the protein sequence. We conjecture here 
that the presence of these "non-sensical" states might act 
as trap states, making a smooth transition towards the 
desired final ground state more difficult and therefore in- 
creasing the time required for an adiabatic evolution. In 
other fields of computer science, physics and chemistry, 
one might also use classical methods or a mean field ap- 
proach to find approximate solutions to be employed as 
educated guesses. For example, in the context of quan- 
tum simulation, a Hartree-Fock solution may be used as 
the initial state for an adiabatic preparation of an exact 
molecular wave function [19.-.22J . 

In quantum mechanics, preparing a desired state is not 
a trivial task [23, [2J| and it is not always possible to de- 
terministically prepare a state from a superposition state 



by measurement. For NP-complete classical problems, 
like the 3-satisfiability problem (3-SAT) [H, [iB] studied 
here, the simplest guess for the initial state is that of 
choosing one of the possible assignments from the so- 
lution space. We present a strategy which allows to de- 
sign experimentally-realizable initial Hamiltonians whose 
ground state corresponds to the desired guess as the ini- 
tial quantum state. This is essential for studying the im- 
portance of choosing a guess state instead of having full 
superposition of states as in CAQC [l|. In an algorithm 
like the one proposed here, both the initial guess and the 
final solution will be states of the computational basis. 
Therefore, the overlap between them is zero, unless one 
has guessed the right solution or the guess is included in 
the subspace of solutions, which is very unlikely. Regard- 
less of this counterintuitive choice, {ip guess {^solution) = 0, 
we show that using this kind of heuristics can be of ad- 
vantage. Even in the case of choosing the initial state 
by random guessing there is potential for outperforming 
CAQC. 

b. Restarting the evolution. Notice that the method 
described here not only allows to begin the evolution from 
a guess state, but also allows for the possibility of restart- 
ing a failed AQC calculation from the measured state. 
This state can be used as a refined guess to restart the 
adiabatic evolution. An adiabatic processes is an ide- 
alized concept because real experiments have to be run 
in finite time and therefore there will always be prob- 
ability of measuring a non-desired excited state which 
does not encode the solution to the computational task. 
The possibility of restarting the quantum evolution using 
the measured state as a guess is a feature which is not 
available in any of the AQC proposals to date, to our 
knowledge. 

The incorporation of heuristics in AQC essentially in- 
volves two modifications to CAQC. We address both 
changes; the first modification involves the design of ini- 
tial Hamiltonians for arbitrary guess states. The second 
modifications involves the change of the time profile of 
the driving Hamiltonian from the linear ramp to a non- 
linear time dependence with a "sombrero-like" time pro- 
ffie (see Fig.[T]). Note that this second change is not the 
main point of thepaper since non-linear paths had been 
proposed before Q. Also, it is not our purpose to ex- 
plore what is the optimal selection for the driving term. 
It must be emphasized that the "sombrero-like" time pro- 
file is an essential feature needed if one is interested in 
the kind of heuristics we describe here, but this is not the 
case in the conventional way of doing AQC where it was 
used for auxiliary Hamiltonian terms ^] . Because of this 
distinctive feature and with the purpose of differentiating 
our heuristic strategy proposed with CAQC, we will refer 
to our method as Sombrero Adiabatic Quantum Compu- 
tation (SAQC). The name should be associated with the 
algorithmic strategy (selection of initial guess, design of 
initial Hamiltonian and sombrero- like profile for the driv- 
ing Hamiltonian) which aims to incorporate heuristics in 
AQC, not only to the use of non-linear paths in AQC. 



The paper is divided as follows: in Section HIl we re- 
view the CAQC approach. Section Hill introduces the ba- 
sic elements of the new implementation, SAQC. Finally, 
in Section HVl we present numerical calculations compar- 
ing the performance of both the CAQC and the SAQC 
algorithms based on the minimum gap, gmin, of their re- 
spective time-dependent Hamiltonians driving their cor- 
responding time evolutions. 



II. CONVENTIONAL ADIABATIC QUANTUM 
COMPUTATION (CAQC) 



The goal of AQC algorithms is that of transforming 
an initial ground state IV'(O)) into a final ground state 
|?/'(t)), which encodes the answer to the problem. This 
is achieved by evolving the corresponding physical sys- 
tem according to the Schrodinger equation with a time- 
dependent Hamiltonian H{t). The AQC algori thm re- 
lies on the quantum adiabatic theorem [25l - l35l |. which 
states that if the quantum evolution is initialized with the 
ground state of the initial Hamiltonian, the time prop- 
agation of this quantum state will remain very close to 
the instantaneous ground state \i/jg{t)) for all t e [0, r], 
whenever H{t) varies slowly throughout the propagation 
time t G [0,t]. This holds under the assumption that 
the ground state manifold does not cross the energy lev- 
els which lead to excited states of the final Hamiltonian. 
Here, we denote by ground state manifold the first m 
curves associated with the lowest eigenvalue of the time- 
dependent Hamiltonian for t e [0,t], where m is the 
degeneracy of the final Hamiltonian ground state. An 
example of to = 2, is shown in Fig. 5 of Rcf. Q- 

Conventionally the adiabatic evolution path is the lin- 
ear sweep of s € [0, 1], where s = t/r: 



H{.s) = (1 - .s)H, 



transverse 



sHf 



(3) 



III. SOMBRERO ADIABATIC QUANTUM 
COMPUTATION (SAQC) 

For SAQC, the time-dependent Hamiltonian can be 
written as: 



H 



sombrero 



(1 - s)H, + hat(s)7Jdr.. 



sH 



/■ 



(4) 



We want the non-degenerate ground state of the ini- 
tial Hamiltonian Hi to encode a guess to the solution, 
and the driving term, H driving, to couple the states in 
the computational basis. The function hat(s) is zero at 
the beginning and end of the adiabatic path; therefore 
Hdriving acts Only in the range s S (0, 1) in a "sombrero- 
like" time dependence (see Fig.[T|), which allows Hi (Hf) 
to be fully turned on at the beginning (end) of the com- 
putation. 



A. Design of the initial Hamiltonian for the guess 
state 

As preparing an arbitrary initial non-degenerate 
ground state for adiabatic evolution is not a trivial task, 
we focus on easy to prepare initial guesses that consist of 
one of the states in the computational basis. The strat- 
egy proposed builds initial Hamiltonians such that the 
initial guess corresponds to the non-degenerate ground 
state of the initial Hamiltonian, as it is required by AQC. 
Additionally, this ground state would be non-degenerate. 

Let us denote the states of the computational basis of 
an N qubit system as \qN) \qN-i) ■ ■ ■ \qi) = \qN ■ ■ ■ qi) 
where g„ G {0,1}. The proposed initial Hamiltonian, 
whose ground state corresponds to an arbitrary initial 
guess state of the form \xn ■ ■ ■ xi), can be written as 



H^ 



N 

E 

n=l 



N 



[Xnl + qn{^-2Xn)] ^'^K„, (5) 

where each Xn is a boolean variable, Xn G {0, 1}, while 
q = ^{I — a^) is a quantum operator acting on the n-th 
qubit of the multipartite Hilbert space Hn (^ ^-N-i (Xi 
• • • <8) "Hn • • • (81 Hi . The operator q^ is given by 



Htransverse (scc Eq. [3] bclow) is usually choscn such 
that its ground state is a uniform superposition of 
all possible 2" computational basis vectors, for the 
case of an n— qubit system. Here, we choose the 
spin states {\qi = 0) ,\q= 1)}, which are the eigenvec- 
tors of (jf with eigenvalues -1-1 and -1, respectively, 
as the basis vectors. Then the initial ground state is 
iV's(O)) = ^ Eg,G{o,i} l9") l^n-i) • • • I92) |gi)- Such an 
initial ground state is usually assumed to be easy to pre- 
pare, for example, by imposing a global transverse field. 
Since each state encodes a possible solution, this initial 
state assigns equal probability to all possible solutions to 
the computational problem. 



qn^ In "S) In~i ^ ■ ■ ■ <S) {q)n ® 



/i 



(6) 



where q is placed in the nth position and the identity 
operators act on the rest of the Hilbert space. 

The states constituting the computational basis, |0) 
and |1), are eigenvectors of CTz with eigenvalues -1-1 and 
— 1, and therefore they are also eigenstates of the oper- 
ator q with eigenvalues and 1 respectively. The logic 
behind the initial Hamiltonian in Eq. [5] then is clear: if 
Xn = 0, then hx,^=o = qn but in the case of Xn = 1, then 

hx„=i = I ~ (jn- 

As an example, suppose one has a four qubit system, 
and one wishes to initialize the adiabatic computation 




0.8 1.0 



FIG. 1: (Color online) The main idea behind Sombrero Adiabatic Quantum Computation (SAQC) is to introduce heuristics 
in AQC, and having the possibility of restarting a failed AQC run from the measured excited state. In order to prepare an 
arbitrary state from any of the 2 possible basis states from the computational basis of the N qubit system, we propose an initial 
Hamiltonian, Hi (see Eq. [5]), in such a way that the desired initial guess state is the non-degenerate ground state of the designed 
initial Hamiltonian. The initial Hamiltonian is diagonal in the computational basis, and so is the final Hamiltonian for the case 
of classical problems such as the NP-complete problems, e.g., random 3-SAT. Since both, the initial and final Hamiltonians are 
diagonal, connecting them via a linear ramp as is usually done in CAQC (see left panel) will not lead the quantum evolution 
towards finding the ground state of the final Hamiltonian. To maintain the initial Hamiltonian uniquely and fully turned on 
at the beginning, f = 0, and the final Hamiltonian uniquely and fully turned on at the end of the computation, t = r, we 
introduce a driving Hamiltonian whose time profile intensity has a "sombrero-like" shape (see right panel) is such a way that 
it only acts during < t < r. Two examples of functions with this functional form are presented, where hati(s) — sin^(7rs) 
and hat2(s) = s(l — s). A desired feature of our algorithmic strategy is the possibility of introducing heuristics, and not that 
of introducing non- linear paths. The latter has been proposed previous publications [3, l3a . |37| . but here is employed as a 
consequence of the algorithmic strategy. 



with the state \x4 = 1,2;3 = 0,a;2 = l,a;i = 0) = |1010) 
which one may choose either randomly or as an educated 
guess to the solution. According to Eq. \E\ the initial 
Hamiltonian for the |1010) guess state should be con- 
structed as 



Hi 1^X4 ' '^xs ' ^X2 ' '^xi 



{i -qd + Qs + U-92) + 9i, 



(7) 



and, clearly 



H, 11010) = ((/ - 94) + 93 + (/ - <?2) + qi) 11010) 
= 0|1010). 



(8) 



In general, the 2^ states of the computational ba- 
sis are all eigenstates of Hi, and it can be easily veri- 
fied that the spectrum of Hi are energies contained in 
{0, • • • , N}. As required, the ground state is also non- 
degenerate. The other states will have an eigenenergy 
which equals their Hamming distance to the ground state 
of the initial Hamiltonian. 



B. Driving Hamiltonian 

The encoding of an educated or a random guess into 
Hi (Eq. [5]) makes both Hi and Hf (Eq. U]) diagonal in 
the computational basis. Therefore, connecting Hi and 
Hf with a linear ramp (as in Eq. [3]), namely, omitting the 

operator Hdriving in the quantum evolution, would yield 
zero probability of obtaining the state that encodes the 
unknown solution to the problem starting from the initial 
guess state. To avoid such a situation, Hdriving must in- 
troduce non-diagonal terms in H sombrero (see Eq. ^ that 
allows the initial state to transform from any arbitrary 
guess into the solution. 

In order to make a fair comparison between CAQC and 
SAQC (see Eq.[3]and Eq.©, we set 



H 



drivi 



H 



transverse 



N 

5^ 



n=\ 



(9) 



in Eq. 21 where q!^ stands for the quantum operator q^ 
acting on the nth qubit of the multipartite Hilbert space 
'Kn ® 'Hat-i ® ■ ■ ■ ® 'Hn ® ■ ■ ■ ^Hi- The operator g^ is 
given by In ® In-i <8) • • • ® {q'^)n (8) ■ ■ ■ <8) /i, where the 
operator q^ = ^{I — a^) has been placed in the nth posi- 
tion, and the /^'s are identity operators. From a physical 
point of view, the Hamiltonians Hdriving and Htransverse 
can be related to a transverse magnetic field. The in- 
tensity of these Hamiltonians is tuned by varying the 6 
parameter. If we set S to be the same for both adiabatic 
algorithms, all the dependence of the transverse field in- 
tensity lies on functions (1 — s) in Eq. |3] and hat(s) in 
Eq. m A reasonable requirement for a fair comparison 



between CAQC and SAQC is that they both provide the 
same average intensity of the transverse magnetic field 
in s e [0, 1]. A choice of hat(s) with the same average 
/(, ha.t{s)ds = J„ (1 — s)ds — 1/2, is hat(s) = 3s(l — s). 

Even though nonlinear evolutions have been proposed 
in previous articles 8, 36, 37], our hat(s) function can 
be as simple or as complicated as desired, as long as 
hat(O) = hat(l) = is fulfilled. There is plenty of room 
to optimize the performance by choosing a more conve- 
nient hat(s) for the adiabatic evolution, we emphasize 
that the additional advantage of SAQC is the possibil- 
ity of choosing an initial guess. In the next section we 
present some results obtained based on one of the simple 
nonlinear function hat(s) — 3s(l — s) and discuss the per- 
formance of both CAQC and SAQC for random 3-SAT 
instances. 



IV. HAMILTONIANS FOR 3-SAT, NUMERICAL 
CALCULATIONS AND DISCUSSION 



In order to provide a proof of concept for SAQC and 
to test the potential usefulness of both random and ed- 
ucated guesses in adiabatic evolution, we performed a 
numerical study on hard-to-satisfy 6- and 7-variable in- 
stances of the 3-SAT problem and compared our results 
with the CAQC approach. Let us now provide a succinct 
introduction to the 3-SAT problem as well as to briefly 
discuss its relevance in the fields of theoretical and ap- 
plied computer science. 



A. Construction of final Hamiltonians for 

satisfiability problems and design of numerical 

calculations 



The K-SAT Problem. Let A 

{ei, 62, . . . , e„, ei, 62, . . . , e„} be a set of Boolean 
variables E = {ei} and their negations E — {e^}. Let 
us now construct a logical proposition P, defined as 

P = AJ(V,tia.)] = KC^, where a, G A, i.e. P 
is a conjunction of clauses Ci over the set A, where 
each clause consists of the disjunction of k literals. 
Proposition P is a K-SAT instance and the solution of 
the K-SAT problem, for instance P, consists of finding 
a set of values for those binary variables upon which 
P has been built (i.e. a bitstring), so that replacement 
of such binary variables for their corresponding binary 
values makes P = 1, namely, proposition P is satisfied. 
3-SAT is a particular case of K-SAT for K=3. 

For example, let us examine the following instance of 
the 3-SAT problem. Let E = {xi,X2,X3,X4,X5,xe} be a 
set of binary variables, and therefore the set of literals is 
A = EU E = {cci, a:2, . • . , Xq,xi,X2, . ■ . , xq}. Consider a 
3-SAT instance specified by the proposition. 



p = 



{xi V 2:4 

(X3 V Xi 

{xi V £2 
[xz V X5 

(X2 V 2:5 

(a;2 V a;3 
(x3 V af4 
(a;2 V X5 
(xs V af5 



VX5) 
Vxe) 
VaTe) 
Vxe) 
Vxe) 
Vxe) 
Vxe) 



A(a;2 

A (X4 
A(X2 

A(afi 

A(X2 

A (afi 
A (:C4 
A(a;3 

A (X4 



VX3 

V:e3 
V:f2 
Vx-3 
Vx-2 
V;r5 
Va;5 

VX5 



V Xi 

V Xq 
\J Xq 

V Xi 

Va;6^ 

V Xq 

V xe" 



A (xi 
A {xi 
A (xi 

A (X2 
A (X2 

A (xi 

A {X2 

A (xi 

A (xi 



VX2 

Var3 
Vx-2 

VX3 

Vx-3 
V Xi 

VX3 
VX3 
M X2 



VX5)A 
VX5)A 
VX6)A 

V X4)A 

V X4)A 
VX5)A 

V XQ)h 

V X6)A 
VX3) 



As this example suggests, finding solutions of even a 
modest 3-SAT instance can become difficult quite easily 
(in this case, P has only one solution: xi = f,X2 — 
1, 2:3 = 0, 2:4 = 1, 2:5 = 0, Xg = 0.) 

3-SAT is an NP-complete problem [ij, [l3: as opposed 
to 2-SAT which can be efficiently solved using a clas- 
sical computer. Consequently, studying the properties 
of 3-SAT is an important area of research, not only be- 
cause a polynomial-time solution to 3-SAT would imply 
P — NP, but also because 3-SAT (due to its polynomial 
equivalence with K-SAT) may be used to model prob- 
lems and procedures in theoretical computer science [38| 
as well as in several areas of applied coinputer science 
and engineering like artificial intelligence [3^, |40| • 

For the purpose of simplifying our discussion, and 
without loss of generality, we randomly generated 3- 
SAT instances with a unique satisfying assignment (USA) 
and their number of clauses to number of variables ratio 
a K, 4.26. This value of a corresponds to the phase tran- 
sition region where hard-to-satisfy instances are expected 
to be found [4l|, |42] ■ For completeness and to avoid any 
kind of bias in selecting this pool of instances, we selected 
2" different instances for every n variable case studied. 
More precisely, to exhaustively study the impact of dif- 
ferent initial guesses with respect to unique solutions in 
the behavior of SAQC, we considered all 64 possible ini- 
tial Hamiltonians Hi (using Eq. [5]) for each one of the 64 
randomly generated 6-variable USA instances. Similarly, 
we built 128 initial Hamiltonians for each 7-variable in- 
stance, one per possible initial guess (see Fig. [2]). The 
instances were selected in such a way that their solutions 
had not only a USA, but also that there was no two in- 
stances with the same solution. 

The generation of our USA 3-SAT instances took sev- 
eral steps, being the first one using the SAT instance 
generator developed by [43]. Unfortunately, as this gen- 
erator does not warranty the production of 3-SAT in- 
stances with unique solutions, we took all generated 3- 
SAT instances and determined, by exhaustive bitstring 
substitution, whether such instances were USA or not. 
We iterated this process until we computed all 6-variable 
and 7-variable 3-SAT instances we needed. 

The number of qubits used in our simulations is smaller 
than state-of-the-art calculations for AQC, such as those 
carried out by quantum Monte Carlo [y, |4J|. Here we 
trade off carrying out few large calculations on many 
qubits for carrying out many calculations on fewer qubits. 



We wish to answer the question: What is the the impact 
of the initial guess in the spectral properties of the time- 
dependent Hamiltonian? To answer this, we explored the 
space of initial guesses in an exhaustive manner for the 
case of 6-variable and 7-variable SAT instances. We run 
a total of 81,920 and 327,680 for 6- and 7- variables re- 
spectively (see Fig. [5]). As shown in the same figure, we 
also numerically explored the importance of the strength 
of the transverse external field for the performance of the 
algorithm. 

Final Hamiltonians Hf are instance-dependent, i.e. 
the structure of each final Hamiltonian depends on the 
particular structure (conjunction of clauses) of each 3- 
SAT USA instance. Our final Hamiltonians Hf comply 
with the property that it must encode, in its ground state, 
the solution to the particular 3-SAT USA instance it was 
designed for [l|, |4^ . The design of the final Hamiltonian 
involves an intermediate step, where a classical cost or 
energy function is constructed for the particular instance 
of interest. Once this energy function is expressed in 
terms of binary variables, it can be easily transformed 
into a quantum Hamiltonian by performing the mapping 
indicated by Eq. ^ where each classical binary variable, 
g„, is transformed into a quantum operator, g„. The en- 
ergy function, Hf^ associated with the final Hamiltonian, 
Hf, can be constructed as a sum of other energy func- 
tions, hci which involve only variables associated with 
one clause at a time. 



H 



f 



Y^hc 



(10) 



Each hc'i is designed such that it is equal to 1 if clause Ci 
is unsatisfied and if the clause is satisfied. Notice that 
the functions hd contribute to the count of unsatisfied 
clauses which defines the spectrum of possible values for 
Hf, with Hf — when all clauses are satisfied. 

Formally, suppose A = {xi, X2, ■ ■ ■ ,Xn,xi,X2, ■ ■ ■ ,Xn\ 
is a set of n binary variables and their corresponding 
negations, P is a 3-SAT USA instance given by P = 
/\j Ci, and each Ci is a disjuction of three elements of A, 



i.e. Ci 



V ajs V Uj with aa,af},a^ e A and indices 



a, P, 7 are natural numbers, not necessarily consecutive. 
Finally, let B = Z1Z2 ■ ■ ■ Zn be a set of n bits to be sub- 
stituted in instance P. Then, hc^ is given by 

0, if substitution of P = Z1Z2 . . . z„ 
in Ua^/ ap W a-y makes Ci ~ 1 

1, if substitution of P = Z1Z2 . . . z„ 
in aa\/ ap y a-y makes Ci — 0. 

To construct such a function for any arbitrary clause 
Ci — Aq V a^ V a-y, it is useful to note that the only 
assignment for which Ci = is when 0^ = 0, a^ = 0, 
and a-y = 0. Therefore, the function h^ by construction, 
should be 1 when Ua = 0, ap = 0, and a^ = 0, and 
otherwise. It can be easily checked that 



he, = (1 - aa)(l - ap){l - Qy), 



(11) 



equals 1 when aa = 0, ap = 0, and a^ = and other- 
wise. Recall that each ap represents a literal and there- 
fore it could be representing the negation of a variable. 
One can always use the identity Xi — \ — Xi to eliminate 
any Xi , and obtain both hd and Hf in terms of the Xi . 
Consider for example the construction of the energy 
function he required for clause C = Sa V x^ V S^ , i.e. C 
is a conjunction of three negated binary variables. In this 
case, C is satisfied by all possible 3-bit bitstrings except 
for 111 and, according to Eq. [TI] the energy function 
assumes the form hd = (1 — Xa){l — a;^)(l — x~f) = 

Ct /3 "Y ' 

As a last example consider a clause of the form C — 
Xa y xp y x^, where C is a conjunction of three non- 
negated binary variables taken from set A. It is clear that 
C will be satisfied by all possible 3-bit bitstrings except 
for 000. According to Eq.[TTl the energy function for this 
clause C is given by he = {I — Xa){l — xii){l — x-y) = 

the first equality of the previous equation, one can easily 
check that he — Q for all possible 3-bit combinations 
except for /ic(OOO) = 1, as expected. 

Once we have the final expression for the final classical 
energy function of Eq. [TUl the final Hamiltonians Hf can 
be obtained using the mapping of Eq.[6l which relates the 
classical binary variables with quantum operators. Since 
we selected only USA instances for our study, each Hf 
has a non-degenerate ground state encoding the unique 
solution of one of our 3-SAT instances with corresponding 
ground eigenvalue equal to zero. The final Hamiltonians 
are the same for both strategies, CAQC and SAQC. 

The numerical results on the dependence of the 
minimum-gap value, gmim as a function of S are shown 
in Fig. [21 Curves were computed by taking the median 
of all bit strings that fulfilled the criteria specified in 
the legend boxes; namely either to produce i unsatisfied 
clauses (UC=i) when substituting the initial guess bit 
string in its corresponding instance, or to be j bit flips 
away from the solution (BF=j). BE represents the well 
known Hamming distance between the solution and the 
initial guess state. We focused on UC and BE because, 
in principle, the notion of closeness of an initial guess to 
the actual solution may be defined with either parame- 
ter. Data corresponding to a fixed value of (5 is a statis- 
tical representation (median) of typical gmin values that 
would be expected for hard 3-SAT instances if the guess 
state belonged to a definite number of UCs or BEs un- 
der an experimental setup using SAQC. Such curves are 
compared with the minimum gap expected for CAQC. 



B. Effects of the variation in the transverse field 
intensity on the minimum energy gap 

The dependence of gmin values as a function of the 
transverse field intensity S leaves open some important 
questions regarding the efficiency of adiabatic quantum 
algorithms, whether CAQC or SAQC. Eor example, what 



is the optimum value of S which minimizes the running 
time of an adiabatic algorithm? How transferable is this 
optimum S value among computational problems? Al- 
though we do not intend to do a thorough study of this 
question in this paper, we would like to give some insight 
into this question and provide a qualitative discussion of 
what kind of results might be expected. 

Eollowing closely the notation from Earhi et al fll] , con- 
sider H{t) ~ H{t/T) ~ H{s), with instantaneous values 
oi H{s) defined by 



with 



H{s)\Ei{s))^Ei{s)\Ei{s)) (12) 



Ea{s)<Ei{s)<---<E2N_,{s) (13) 



where 2^ is the dimension of the Hilbcrt space, and N the 
number of qubits or equivalently the number of binary 
variables in the SAT instance. According to the adiabatic 
theorem, if the gap between the two lowest levels, Ei{s)~ 
Eq{s), is greater than zero for all < s < 1, and taking. 



r > 



£ 



5r, 



with the minimum gap, gmin, defined by. 



gmzn = niin {Ei{s) - Eo{s)), 

0<s<l 



and £ given by. 



.dH. 



f = max|(i?i(s)| — |i?o(s))|, 
o<s<i as 

then we can make the normed overlap 

|(i?o(. = l)|^(r))| 



(14) 



(15) 



(16) 



(17) 



arbitrarily close to 1. In other words, the existence of 
a nonzero gap guarantees that \ip{t)) remains very close 
to the ground state of H{t) for all < t < r, if r is 
suSiciently large. 

Even though we are aware of the new and more strin- 
gent conditions for adiabaticity |26l - l35| and that Eq.lMlis 
just one of the inequalities to guarantee adiabatic evolu- 
tion (though there is still lack of a sufficient and necessary 
condition according to Ref . [46| ) , we will base our discus- 
sion on Eq. [14] to illustrate that there is nothing anoma- 
lous in employing the additional Hamiltonian term in the 
full time-dependent Hamiltonian for SAQC. As well as in 
CAQC, the algorithmic complexity relies again in avoid- 
ing an exponentially narrowing of gmin- Along the way 
we find an important observation about the scaling of the 
running time as a function of the parameter S. 

Let us first determine an upper bound for Ssaqc 
in Eq. [141 Consider the Hamiltonian in Eq. [4l with 
hat(s) = 3.s(l — s) since this was the functional form used 
for our numerical calculations. We already discussed, at 
the end of Sec. lIII Al that the spectrum of Hi is contained 
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for a total of (128USA) x (128guesses) x (20values of &) = 327,680 SAQC settings. 



-e-BF=0 (128) 


-)<-BF=l (896) 


-^BF=2 (2688) 


-^BF=3 (4480) 


t(^BF=4 (4480) 


BF=5 (2688) 


^BF=6 (896) 


-^BF=7 (128) 


Conv. AQC (128) 




2 3 4 5 5 7 

8 (intensity of transverse field) 
(a) 



0.9 

0.8 

0.7 

I 0.6 

I 0.5 

E 
1 °'' 

0.3 

0.2 

0.1 

0, 



-a-UC=0(128) 
-^UC=1 (1039) 
-^UC=2(2661) 
"' UC=3(3786) 
^(eUC=4(3734) 
UC=5 (2685) 
-|>-UC=6(1458) 



UC=7 (609) 
UC=8 (222) 
UC=9 (50) 
UC=10(9) 
UC=11 (3) 
Conv. AQC (128) 




2 3 4 5 6 7 

8 (intensity of transverse field) 
(b) 



FIG. 3: (Color online) Summary of the 327,680 calculations for 7 variable 3-SAT instances of minimum-gap median values 
as a function of the transverse field intensity within groups sorted by (a) number of bit flips, BF, representing the Hamming 
distance between the initial guess and the solution and (b) number of unsatisfied clauses (UC). Plots include the CAQC (see 
Eq. [3| results for the same 128 different unique-satisfying assignment 3-SAT instances randomly generated for the case of 7 
variables. Values in parentheses correspond to numbers of data points which contributed to the value of the median plotted in 
each curve. Results for 6 variable instances (not shown) are similar to the ones shown for 7 variable instances. 



in {0,1,--- ,N} and, similarly, it can be easily shown 
that the spectrum of Hdriving (see Eq. ^ is contained 
in {0, 1, • • • , SN}. On the other hand, the spectrum of 
the final Hamiltonian, Hf, is instance dependent, and 
its construction guarantees that the maximum eigenvalue 
would be M which denotes the total number of clauses. 
This eigenvalue M would only appear in case we had an 
assignment which violates all of the clauses. Using these 
spectra upper bounds, we can establish an upper bound 
for SsAQC in Eq. [H i.e. 



SsAQC = max I {El (s) \-^\Eo (s)) | 

0<S<1' US 

= max I {El {s)\Hf - H, + iSH driving 

- &6sHdriving\EQ{s))\ 

< m&^{\{Ei{s)\Hf\E^{s))\ + \{Ei{s)\H,\Eq{s))\ 

+ 'i\5{Ei{s)\Hdriving\Ea{s))\ 
+ &\5{Ei{s)\Hd„ving\Eo{s))\) 

<M + N + i\5\N + Q\5\N = N{a + 1 + 9\S\) 

(18) 

Where we have used the triangle and Schwartz inequality 
and also the fact that M = aN, with a close to 4.26 in 
this particular study. We can see that in the worst case 
scenario, Ssaqc scales linearly with the number of vari- 
ables N, and linearly with the intensity of the magnetic 
field, SsAQC — 0{\S\N). A similar analysis gives also 
that ScAQC — 0{\S\N), and therefore, we showed that 
for SAQC, not surprisingly, the algorithmic complexity 
also relies on the scaling of gmin ■ 

An interesting observation arise by analyzing the lin- 
ear scaling of £ with 6 and using the numerical results for 
the dependence of the typical minimum gap values as a 
function of S. There seem to be at least two distinguish- 
able regimes for the dependence of gmin on S for both 
CAQC and SAQC (Fig.JS]). For relatively small values of 
6 G [0.5,1.5], gmin scales approximately linearly with S 
and therefore gmin ~ ^^ ■ Since the running time is given 
by Eq. ll4[ and £ ^ S, the running time t decays inversely 
proportional to 6 within this linear regime. However, for 
large values of S, in the 'stationary' regime where gmin 
is almost constant, increasing field intensity through 5 
would make both algorithms less efficient as running time 
T would increase roughly linearly with S. 

Both CAQC and SAQC would benefit from an increase 
in the transverse field for small values of 6, but notice that 
gmin values for SAQC are more sensitive to S, and soon 
become better on average than those for CAQC (Fig. [3|) . 
According to the previous discussion about running time 
as a function of S, it would be ideal to choose i5 near the 
end of the linear regime; in our calculations, S somewhere 
in the interval (1,2). Further studies concerning the op- 
timum value of (5 as a function of the number of binary 
variables are needed, but we chose 5 = 1.5 for our anal- 



ysis on the performance in SAQC and CAQC described 
in the following section. 



C. Performance comparison betw^een SAQC and 
CAQC 

The data sorted with respect to BFs and UCs shows 
an increase of the minimum gap, gmin, as the Hamming 
distance from the initial guess to the solution decreases; 
the trend for UC is less apparent (see Fig. ^. Comput- 
ing the number of UCs produced by a given initial guess 
can be done in polynomial time on a classical computer. 
Unfortunately there is no way to determine a priori how 
many bit fiips the guess is from the solution, as that 
requires knowledge of the solution itself. Additionally, 
Fig. \M,^) shows that the SAQC implementation, using 
Hdriving ^s defined in Eq. 1^1 does not necessarily favor 
states with low values of UC, but rather gives a homoge- 
nously distributed success probability between 25-45%, 
for 5 — 1.5. This is in accordance to the observation that 
solving 3-SAT hard instances is not necessarily guided by 
minimizing the number of UCs [47[ . Given the above sce- 
nario, we analyzed the likelihood of better performance 
by choosing initial guesses at random. 

In the following discussion, we use the term signifi- 
cantly better initial guess to mean an initial condition 
that leads a SAQC algorithm to be at least twice as fast 
as CAQC, i.e. running times for CAQC, tcaqc, and 
SAQC, TsAQC, are such that tcaqc > 2 tsaqc or, equiv- 

alently, g^^^'^ > V^gZfn'^: assuming fcAQC = ^SAQC- 
For S — 1.5, choosing an initial state at random yields a 
probability greater than 50% of having g^j„ > gmin 
(squares) as shown in Fig. IDJb). Moreover, the proba- 
bility of significantly better performance, i.e. tcaqc ^ 
'^ TSAQC is ~ 35% (triangles). With the intention of pre- 
dicting the performance of the SAQC protocol in the 
limit of large n, the third curve (circles) was produced us- 
ing the following rationale: for USA instances, the num- 
ber of bit configurations with a given value of BF — m 
follows a binomial distribution (J^) . In the limit of large 
n, the likelihood of choosing a state in the central region 
of the binomial distribution is the highest. This obser- 
vation led us to concentrate on the performance of the 
most populated instance subsets, those that correspond 
to BF = 3,4 for 7 variables. Here, the probability of 
significantly better performance is close to 40%. 

Finally, we propose an algorithm based on SAQC. An 
initial guess is chosen either at random or by applying 
expert-domain knowledge and then encoded into the ini- 
tial ground state of Hi (Eq. [5]). An adiabatic passage 
based on SAQC is then performed either in serial or in 
parallel, depending on the availability of quantum hard- 
ware resources (see Fig.[Sl). As an example of the poten- 
tial usefulness of our algorithm, recall from Fig. 0] that 
the probability of significantly better performance using 
SAQC is 39% for 6 = 1.5. 

One way to employ the probabilities we obtained from 
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FIG. 4: (Color online) (a) P{gmin — ^^9min l^^ ~ ^) i^ ^^^ conditional probability of choosing a state with UC= n and 
a SAQC minimum-gap large enough so that the performance of the SAQC is significantly better (at least twice faster) than 
the CAQC. The results were obtained for S — 1.5. Panel (a) shows that there is no correlation between the number of violated 
clauses and the Qmin of the SAQC algorithm, for the hard-to-satisfy instances randomly chosen for this numerical study. Panel 
(b) shows the probability of choosing an initial state at random and satisfying the condition specified in the legend, for different 
values of the transverse field intensity, 0.5 < S < 10. The conditional probability P{g^^'^ > V2gZf^'^\BF = 3,4) (triangles) 
aims to predict the performance of the SAQC algorithm in the case of large number n of qubits. In this limit, an initial state 
chosen at random will have with high probability a Hamming distance BF ~ n/2, given that they are binomially distributed, 
i.e., the number of n bit strings with BF = m is equal to (") 



n bits 




Computational 
basis 




Start over again 



Pick a different 



gues^ 





SAQC 




Guess 1 ^ 












o 

CD 

O 






FIG. 5: (Color online) Implementation of an SAQC algorithm either in parallel or in serial. The algorithm begins by choosing 
a state from the computational basis. For each chosen initial state, an initial Hamiltonian is prepared according to Sec. IIII Al 
Next, choose an ideal time assuming a CAQC protocol will be run, and use that as a reference to run the SAQC protocol twice 
as fast. If only one AQC computer is available, one can still use the probabilistic speed up obtained in SAQC to run for example 
two adiabatic protocols instead of one, in serial mode. Once the first SAQC calculation is finished, one can efficiently check 
whether or not the result is a solution. In case that it is not a solution, one can submit an additional calculation, either randomly 
selecting another initial guess state or using the measured excited state. We call this latter option, "quantum heuristics" since 
the outcome of the near-adiabatic quantum evolution is used to refine the initial guess state for further experiments. In the 
case of having several adiabatic quantum computers at hand, one can do the same initial procedure of selecting guesses, but 
now submitting a different guess to a different node and running on each node twice as fast. 
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our numerical simulations in a more concrete scenario is: 
suppose one is assigned the task of using an adiabatic 
quantum computer and assume one uses Eq. [T4|or any of 
the more stringent conditions for adiabaticity [26l435l | to 
estimate for how long one may need to run an algorithm 
under the CAQC paradigm. Moreover, suppose that this 
estimated running time required to remain in the ground 
state with a high success probability is tcaqc = 2 days. 
Using the numerical results presented in this paper, one 
can opt for performing the same task using SAQC as 
follows; suppose we have absolutely no information about 
the problem as it was in the case of this numerical study 
with random 3-SAT USA instances. Instead of running 
the CAQC algorithm for 2 days, pick a guess state at 
random from all the possible assignments and then use 
it as initial state for SAQC and run the algorithm with 
TSAQC = 1 day. According to Fig. |4l the probability 
of having picked a state whose performance is as good 
as CAQC is 39%, for S — 1.5. If after measurement at 
the end of the first day the result is not a solution, we 
still can pick another state at random and let it run for 
one more day. By now the probability of having picked a 
state with the same performance as the CAQC in the two 
trials equals 63%. Note that in this simple probability 
calculations we are not taking into account the fact that 
even in the case where the state selected for the first run 
was not one of the 'ideal' ones (states corresponding to 
the 39 % of guesses for the results presented in Fig. |4]), 
we still have a very good chance that the 'non-ideal' state 
still delivers a right answer after the first measurement. 
This probabability will depend of course on how close is 
the chosen state from the set of 'ideal' ones. 

Consequently, the execution of two SAQC algorithms 
in serial would take at most as much time as the exe- 
cution of only one CAQC algorithm. By allowing us to 
choose two guesses to run in the same time as one case 
in CAQC, the probability of choosing a significantly bet- 
ter initial guess in these two SAQC executions increases 
from 39% to 63%. Furthermore, even when no signifi- 
cantly better initial guess is chosen and the process is 
not guaranteed to be fully adiabatic, there is still some 
probability that we measure the correct solution at the 
end of both executions. 



V. CONCLUSIONS 

In summary, we propose an algorithmic strategy which 
incorporates heuristics in adiabatic quantum computa- 
tion. In particular, we study one of the most basic heuris- 
tic strategies consisting of initializing the computation 
with a desired initial state chosen by physical intuition, 
an educated guess, classical preprocessing of the prob- 
lem, and/or by randomly choosing one of the possible 
assignments. This method allows to bridge powerful clas- 
sical techniques such as heuristic optimization routines 
and/or mean-field calculations to obtain approximated 
solutions of the problem and use them as initial guess 



states for SAQC. The strategy presented allows for a par- 
allel and/or a concatenated scheme. The parallel setup 
might be helpful if several adiabatic quantum platforms 
are available in which several guesses can be run simul- 
taneously, one guess to run in each one of the adiabatic 
processing units. The idea of the concatenated scheme is 
to restart a failed adiabatic evolution with the measured 
excited state. Assuming a near- adiabatic trajectory, the 
measured excited state can be taken as a refined guess 
from which the quantum evolution is restarted. Notice 
this is not possible in the conventional approach which 
would restart the evolution from the full superposition, 
attributing equal probability to all the states, therefore 
"erasing" the information gained in the previous near- 
adiabatic evolution. Neither of the two features pro- 
posed above are possible in any of the different adiabatic 
quantum computation proposals to date. In addition, all 
the modifications proposed related to different adiabatic 
paths, auxiliary perturbations |48| can also be explored 
in the context of S AQC, w hich is also suitable to study 
quantum problems J49l453j . 

The numerical study performed in this paper is a proof- 
of-principle to explore the importance and consequences 
of starting the adiabatic evolution with a guess state and 
to illustrate that incorporating this kind of heuristics in 
AQC is possible and might be advantageous. Our numer- 
ical simulations show that starting the adiabatic evolu- 
tion with a guess state which has a zero-overlap with the 
solution, {(l)guess\4>soiuUon) = 0, is uot a big concern. On 
the contrary, even when there is no hope to make an ed- 
ucated guess and selection of the initial state at random 
is the only available alternative, we obtain that approx- 
imately 40% of states might allow running the quantum 
algorithm at least twice as fast when compared to CAQC. 
This possibility of running the algorithm for shorter times 
but with several trials brings also additional advantages 
of getting the right answer in any intermediate measure. 
Moreover, these shorter runs would be less affected by 
decoherence effects. 

Even though the procedure used for the performance 
comparison, CAQC vs. SAQC, is "reasonable" since it is 
based on the widely-used minimum gap criteria and its 
connection with the algorithm run-time (see E q . [T4l) , we 
are aware of the limitations of this analysis |26l - t35l l . We 
want to stress that the present numerical results are only 
encouraging indicators that heuristics in AQC might be a 
valuable algorithmic strategy for AQC, given the strong 
dependence of the value of the minimum gap as a func- 
tion of the initial guess chosen. It is not our purpose to 
claim superiority of SAQC but to introduce the approach 
and the motivation behind it. For a more rigorous com- 
parison of both schemes, CAQC and SAQC, we suggest 
numerical experiments which are not meant to be ex- 
haustive but preferably involving larger size instances, 
and to explore different problems other than random 3- 
SAT. For example, we are interested in performing these 
studies for relevant instances of our recently developed 
AQC proposal for protein folding [7| and for adiabatic 
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preparation of molecular ground states [19| where we ex- 
pect the mean-field (Hartree-Fock) solution to be a better 
guess than a full superposition. Numerical propagation 
of the time-dependent Schrodinger equation, instead of 
inspection of the minimum gap after the Hamiltonian 
diagonalization, for these cases will provide a realistic 
simulation of the quantum computation. 

Open questions to be explored further is the connec- 
tions between SAQC, quantum phase transitions [54 
entanglement (Sg and the effect of local minima 
The performance of the adiabatic algorithm in the limit 
of large n is still an open question I4J) [531 which needs 
to be explored in the context of SAQC. It is not obvious 
that the same observations and conclusions/observations 
mentioned above for CAQC will hold for SAQC as well. 
For example, we think that SAQC might be a stra tegy 
to avoid the local minima traps described in Ref. |13| . 
From the study of the spectral properties of the ran- 
dom 3-SAT instances studied, we found that some initial 
states have a considerably larger gap and while others 
show a considerably smaller gap when compared with 
CAQC. Since the initial state in CAQC is a full super- 
position of all the possible states or solutions, including 
the ones with a large gap and others with small gap, 
we conjecture here that having a full superposition will 
not necessarily be the best choice, given that the pres- 
ence of the states with small gaps could slow down the 



quantum evolution. An analysis beyond the gap criteria 
would be needed to test this conjecture. Solving the time- 
dependent Schrodinger equation for the entire evolution 
is the most straightforward, yet numerically-challenging 
approach. In SAQC, the probability of obtaining these 
trap states can be avoided and even in the case of a failed 
evolution, the concatenated scheme may help to restart 
using a better guess. In contrast, in CAQC, the full su- 
perposition including states with considerably small gap 
might result in a bottleneck for the dynamics towards a 
successful computation. Further studies need to be done 
to verify this is indeed the case. 
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